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We consider the dimerized spin-i Heisenberg chain with 
spin hexameric distortion of the exchange pattern and 
study the zero-temperature phase diagram in the pa- 
rameter space ( Ji, J2, J3) by continuum-limit bosoniza- 
tion approach and the exact diagonalization method. The 
phase diagram is rich and has two gaped dimer phases. 



We obtain an estimate of the critical line separating the 
different gapped dimer phases by the bosonization ap- 
proach. The existence of the transition line and the dif- 
ference between dimer phases is checked numerically. 
The behavior of the energy gap and the dimer order pa- 
rameter supports the exact location of the gapless line. 
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1 Introduction The critical gapless phases emerging 
on the border lines separating two gapped phases in low- 
dimensional spin systems has been the subject of studies 
for decades. The spin S = 1/2 Heisenberg chain with 
dimerization and frustration, is a primer and well studied 
example of a such system which shows a gapless phase in 
the ground state phase diagram II]|2][3]|4l . After the semi- 
nal paper by Martin-Delgado et.al [5i], great attention has 
been focused on the studies of the same phenomena in the 
other wide class of low-dimensional magnets such as the 
spin ladders. ^l7ll8ll^flWTTfnifT3lfT4l Recently, the gap- 
less phases on the border of different massive phases has 
been discussed in two-dimensional spin systems 11511161 

no. 

In this paper, we address this problem within the 
slightly different framework, namely we study the stabil- 
ity of the explicitly dimerized (gapped) Heisenberg chain 
towards the perturbation caused by the different commen- 
surate modulation of the exchange pattern which higher 
period itself causes a gapped state however breaks the 
favored by the explicit dimer order. 

The Hamiltonian of the model that under consideration 
is defined as 



where 



N 



ff = J^(l-(-l)' , <5)S n .S n+1 , 



(2) 



is the Hamiltonian of the dimerized chain with strong odd 
links (5 > 0) and 



N 

H x = AJ g{n)S n ■ S n+1 



where 



g(n) = 1 + cos 7m + 2 



cos(-n) 



2tt . 
cos — n) 
v 3 ; 



(3) 



(4) 



H — H + H\ , 



(1) 



is the perturbation Hamiltonian of the dimerized chain 
with commensurate and modulation of the exchange on 
each 6-th link which increase (decrease) at A > (A < 0). 
Here S„ is the spin S = i operator at site n and the chain 
consists of TV = QNo sites. It is straightforward to obtain 
the Hamiltonian Eq. ([TJ, describes a Heisenberg chain with 
the following hexamer modulation of the spin exchange 
(see Fig.Q]). 

J(l) = 7(3) = 7(5) = 7(1 + 6) = 7i 

7(2) = 7(4) = 7(1 -5) = 7 2 , (5) 

7(6) = 7(1 - S + 6A) = 7 3 . 
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Figure 1 (Color online.) Schematic representation of spin 
chains with hexameric modulation of spin exchange which 
is considered in the paper. 



In our recent paper 0181 . we have studied the magnetic 
phase diagram of the model in the limit of strong exchange 
on odd links. It has been shown, that the presence of ad- 
ditional modulation leads to the dynamical generation of 
two new energy scales in the system and to the appearance 
of two additional quantum phase transitions in the ground 
state of the system. These new gaps appear at finite mag- 
netization and show themselves in the presence of two new 
magnetization plateaus at magnetization equal to J and | 
of the saturation value. 

In this paper, we study the effects caused by the ad- 
ditional distortion of the exchange pattern with period six 
on the ground state properties of the dimerized chain at 
zero magnetization. We show that in the case of frustrating 
character of hexamer distortion, with increasing A a quan- 
tum phase transition takes place in the ground state of the 
system. The transition corresponds to the reestablishment 
of the gapless critical phase on the border line, which sep- 
arates two dimerized phases with shifted in respect to each 
other on one lattice unit dimer order. 

The paper is organized as follows: In the next section, 
the bosonization method is applied and the transition line 
is determined analytically. In the section 3, the results of a 
numerical experiment are presented. Finally, we will dis- 
cuss and summarize our results in section 4. 

2 Bosonization In this Section, we use the continuum- 
limit bosonization treatment to study the ground state 
phase diagram of the model (Q3. To obtain the continuum 
version of the Hamiltonian, we use the standard bosoniza- 
tion expression of the spin operators |[T9l 

1 ~ , (-l) r: 



S ,, 



d x 4>(x) 



ira 



■ sin(v27r</>) , 



S% = cos(\/2^6 



1 



(— 11™ 

^ — '— wa.(y/2K6) 
na 

(-1)" 
7ra 



sin(v 27r</>) 



Here <j>(x) and 9(x) are dual bosonic fields, d t <j) 
and satisfy the following commutation relations 

[<j>{x),6{y)}=iG{y-x), 
[<Kx),0(x)]=i/2, 



(6) 
(7) 

(8) 

v s d x 9, 

(9) 



and ao is the constant of the order of lattice unit. Using 
©-(O, one gets the following bosonized expression for the 
alternating part of the "dimerization operator" 



S„ • S n+ i ~ —5 — cos(V2n<t>(x)), 



ir*a 



(10) 



and finally (see for details Ref. |[T8l ) the following 
continuum-limit bosonized version of the Hamiltonian 

H Bos = J dx{^[(d x( f>) 2 + (d x 6) 2 } 
A 



27rao 



2^) }, 



(11) 



where 



A = J{S-\) = -(3J 1 -2J 2 -J 3 ), (12) 
6 

and the velocity of spin excitations 

u = «o J(l + 6A) = ao Jeff- 

From the exact solution of the sine-Gordon model, it is 
known that for arbitrary finite Aq the excitation spectrum 
of the Hamiltonian Eq. ( fTTT i is gapped and consists of soli- 
tons and antisolitons with mass M and soliton-antisoliton 
bound states ("breathers") with the lowest breather mass 
also equal to M (20 ,21 1. The soli ton mass M, which de- 
termines the gap in the excitation spectrum and is con- 
nected with bare model parameters (Aq and J e ff) as fol- 
lows M = J eff C{A /J eff ) 2 / 3 where C ~ 1.4 E2- 

The gapped character of the excitation spectrum results 
to suppression fluctuations in the system and the 4> field is 
condensed in one of its vacua ensuring the minimum of the 
dominating potential energy 



(0) 







at A < 
at A > 



(13) 



The vacuum expectation value of the cosine field in the 
gapped phase in the weak coupling is given by 



(COS ' 



^4>) * {M/ J eff ) 1/2 



(14) 



while in the strong coupling, at \M\ > J it becomes of the 
unit order ll23l . 

In absence of the hexamer distortion, at A = (i.e. 
J 2 = Jz = */(l — 5)), Aq = 5 J > 0, the excitation spec- 
trum is gapped and the field <f)(x) is pinned with vacuum 
expectation value ((f) = 0. Using the bosonized expres- 
sions for spin operators one can easily get the on-site spin 
order is strongly suppressed, while the link-located dimer 
order 



(S„ • S n +i) = Const + (- 



-ir4 



(15) 



shows the long-range ordered for dimer order and the max- 
ima of the dimerization functions located on even links. 
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At 5 = 0, i.e. Ji = J 2 = J, J 3 = J(l + 6A) again 
the excitation spectrum is gapped, but in this case with the 
bare mass Aq = —A J < 0. In this case for A < 0, (f>(x) is 
pinned with vacuum expectation value (<f>) = and for A > 
0, c/>(x) is pinned with vacuum expectation value (<fi) = it. 
As a results, fluctuations of the on-site degrees of freedom 
are fully suppressed , while the link-located dimerization 
function shows a long range order 

(S„ • S„+i) = Const + (-1)" sign{\) (16) 

At A < 0, exchange on each 6-th weak link becomes 
weaker, the dimer order introduced by the hexamer distor- 
tion of the spin exchange is once again given by Eq. il5[ 
and therefore as in the case with S > maxima of the 
dimerization functions are located on even links. 

At A > the dimer order is pinned with the strongest 
links and the minimum of the energy is realized with 
dimer order which shows maxima on even, including the 
strongest each 6-th link. Therefore, in marked contrast 
with A < case, at A > the hexamer distortion plays 
the role of frustration with respect to initial explicit dimer 
order. 

Therefore, at finite S > and A > 0, in the gapped 
phase, the long range ordered dimerization pattern has 
maxima on odd links at Aq > and on even links, at Aq < 
0. In the former case, the distribution of dimers coincides 
with the order of non-disturbed dimerized chain, while in 
the case of So = 0, the order determined with dominant 
strong links. These two gapped phases with shifted on one 
lattice unit with respect to each-other pattern of dimers are 
separated by the critical line where Ao = i.e. 

3Ji - 2 J 2 - J 3 = 6 J(S - A) = , (17) 

and the gapless critical phase which described by the gaus- 
sian free field is realized. 

In order to investigate the detailed behavior of the 
ground state phase diagram and to test the validity of 
the picture obtained from continuum-limit bosonization 
treatment, below in this paper we present results of numer- 
ical calculations using the exact diagonalization for finite 
chains. 

3 Numerical experiment In this section, we study 
the ground state phase diagram of the dimerized Heisen- 
berg chain with spin hexameric distortion of exchange by 
performing a numerical experiment. One of the remark- 
able ways in field of numerical experiments is known as 
the exact diagonalization (ED) method. We have used the 
ED method to diagonalize numerically the Hamiltonian 
([]]) with periodic boundary conditions. To determine the 
ground state phase diagram of the model, we have cal- 
culated the spin gap and the dimer order parameter as a 
function of coupling exchange on odd links, J\, for finite 
chains with lengths up to TV = 24 and different values of 
exchanges J 2 , J3. 

In Fig. [2] we have presented results of our numer- 
ical calculations on the energy gap for the values of 




1/2 3 

Figure 2 (Color online) The energy gap versus J\ for the 
dimerized chains with lengths TV = 12, 18, 24 and different 
values of the exchanges (J 2 , J3) = (a) : (1.0, 1.0), (b) : 
(2.0, 5.0) , (c) : (3.0, 9.0) and (d) : (1.0, 4.0). 



the exchanges parameter corresponding to (J2,Js) = 
(a) : (1.0, 1.0), (6) : (2.0, 5.0), (c) : (3.0,9.0) and 
(d) : (1.0, 4.0) for different chain lengths TV = 12, 18, 24. 
We have recognized the energy gap in finite chains as the 
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difference between the energies of the ground and first ex- 
ited states. It is clearly seen that the spectrum of the model 
is gapped at J\ = 0. As soon as the J\ increases from zero, 
the energy gap decreases and vanishes at the critical value 
Jf . It is obvious that for finite systems, the energy gap is 
always finite and only vanishes at certain value J c { in the 
thermodynamic limit TV — > oo. By further increasing of 
J\, the energy gap reopens and behaves almost linearly in 
respect J\ . Therefore, the peculiar behavior of the energy 
gap shows that the system can be found in two different 
gapped phases by tuning the exchange values. 

At J i = 0, the system reduced into N/2 pair of spins 
in the singlet state. It is completely natural to expect for 
finding the ground state of the system in a phase with 
dimerization on even links, so called Dimer-I. On the other 
hand, in the limit of very strong exchange on odd links, 
J-2 = J3 = 0, the ground state has the structure of sin- 
glets on odd links, so called Dimer-II. So, the competition 
between two terms in the Hamiltonian should be lead to a 
quantum phase transition between these two dimer phases. 
It is known that due to the nature of a dimer phase, trans- 
lation invariance symmetry is broken by one unit cell of 
the lattice 11111 . The structure of the dimer phase can be 
obtained by studying the dimer order parameter defined as 

2 N/2 

D = jjY,< Gs|S 2 „-i.S 2n - S 2 „.S 2n+1 |Gs >, (18) 

n=l 

where the notation (Gs\...\Gs) represents the ground state 
expectation value. It is expected that the dimer order pa- 
rameter, D, takes the values | and ^ in the saturated 
Dimer-I and Dimer-II phases respectively. We have pre- 
sented our numerical results on the dimer order parame- 
ter in Fig.[3]for finite chains with lengths N = 12, 18, 24 
and exchange values (J 2 , J3) = (a) : (1.0, 1.0), (6) : 
(2.0, 5.0), (c) : (3.0, 9.0) and (d) : (1.0, 4.0). As is clearly 
seen, in the absence of Ji, the dimer order parameter is 
positive and the ground state of the system is in the satu- 
rated Dimer-I phase. As soon as the J\ is added and in- 
creased from zero, induced quantum fluctuations reduce 
the dimer order parameter D and it reaches to zero value 
at the critical point Jj\ We found that, taking zero of D at 
the critical value Jf is size-independent which shows its 
validity in the thermodynamic limit N — > 00. By more 
increasing J±, the dimer order parameter becomes negative 
which is the indication of the long-range dimer ordering on 
odd links, namely Dimer-II phase. In this region, the neg- 
ative value of D enhances by increasing the exchange J\ 
and it tends to the saturation value of the Dimer-II phase. 

To find a better physical picture of the gapped dimer 
phases, we have a look into the microscopic behavior of 
the system by using our numerical experiment. Since in a 
unit cell of a hexamer chain (Fig. [T}, there are three differ- 
ent links we have calculated the dimerization on each link 




Figure 3 (Color online) The dimer order parameter versus 
J 1 for the dimerized chains with lengths N = 12, 18, 24 
and different values of the exchanges (J 2 , J3) — (a) : 
(1.0,1.0), (6) : (2.0,5.0) , (c) : (3.0,9.0) and (d) : 
(1.0,4.0). In the inset, the functions /(l),/(2) and /(3) 
are plotted. The dashed line is only a guide for the eye. 
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J./J =2/3+J,/3J. 



Dimer-II 



V J 2 



Figure 4 (Color online) The phase diagram of the dimer- 
ized Heisenberg chain with hexameric distortion of the 
exchange pattern. Solid squares are the continuum-limit 
bosonization data, where solid circles represents the exact 
diagonalization results extrapolated into the infinite chain 
limit. 



defined as 

/(l) = < S6n-5-S6n-4 _ S6ri-4-S6n-3 >) 
/(2) = < S6n-3-S6ri-2 _ S6n-2-S6n-l >j 
/(3) = < S6 n -l.S6„ — S6n-S6n+1 > • 



(19) 



In the Inset of Fig. [3] We have plotted numerical results 
on these microscopic functions versus the exchange J\ 
for chains with lengths N = 12, 18, 24 and exchanges 
(J 2 ,J 3 ) = (a) : (1.0, 1.0), (6) : (2.0, 5.0), (c) : (3.0,9.0) 
and (d) : (1.0,4.0). It is clearly seen that the values of 
functions /(l),/(2) and /(3) take zero at different val- 
ues of J i (dashed lines). By considering the strength of 
exchanges, expected that the first the dimer on the weak- 
est links takes zero and finally, the strongest one becomes 
zero which behavior is completely seen in our numerical 
results. 

We have to mention that for other values of the ex- 
changes, we did our numerical experiment as well and 
found the same ground state magnetic phase diagram con- 
tains of: (I.) gapped Dimer-I phase in the region Ji < 
(II.) gapped Dimer-II phase in the region Ji > Jf (Fig.|4|i. 
In this figure, solid squares are the critical points obtained 
using the continuum-limit bosonization and solid circles 
represents critical points obtained by numerical experi- 
ment. The comparison of numerical results with predic- 
tions from the continuum-limit field theory shows well 
agreement. Besides the above calculations, to obtain the 
type of the mentioned quantum phase transition between 
dimer phases in our model, we have implemented the al- 
gorithm to find the ground state energy. A very important 
indication of the first order phase transition is the discon- 
tinuity in the first derivation of the ground state energy at 



•a 
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0.05 - 



0.00 



(a) J =3 , J 2 =2.5 
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N=24 
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-a 
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0.05 - 



(c) J =3 , J =3.5 




1 2 \ ' ' " 3 4 5 



Figure 5 (Color online) The first derivation of the ground 
state energy as a function of the control parameter J3 
for the dimerized chains with lengths N = 12, 18, 24 
and different values of the exchanges (Ji, J2) = (a) : 
(3.0,2.5), (b) : (3.0,3.0) and (c) : (3.0,3.5). In the in- 
set, the second derivation of the ground state energy are 
plotted. 



Copyright line will be provided by the publisher 



6 



the quantum critical point. Using the numerical Lanczos 
method, we have calculated the ground state energy for 
chain sizes N = 12, 18, 24 and plotted the first derivation 
of the ground state energy as a function of the control pa- 
rameter J3 in the Fig. [5] The results of the first derivation 
show absence of any discontinuity and therefore in agree- 
ment with bosonization studies indicates on the continu- 
ous character of the phase transition. Also, we have calcu- 
lated the second derivation of ground state energy in the 
inset of Fig. [5] It is clearly seen that the height of peak 
enhances by increasing J3 and the second derivation of 
ground state energy will be diverged in the thermodynamic 
limit N — > 00 that is in agreement with the continuous 
phase transition. The shift in the position of peaks is result 
of the finite size effect. 

4 Conclusions In this paper, we have studied the 
ground-state properties of the dimerized Heisenberg chain 
with spin |, which has hexameric distortion of the ex- 
change pattern. We have determined the existence of a 
gapless line in the quantum phase diagram by the bosoniza- 
tion technique in the continuum-limit approach. Moreover, 
the exact diagonalization method has indicated that the 
ground state has two gaped phases, so called the Dimer-I 
and Dimer-II phases. Our numerical results are in well 
agreement with analytical results. 
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\ "Km-TX nmocra V 

\UIil \ XUIIltJHjCLj 


J] 


\Dnii_ \ e taj- 


? 


\Dmi\xij- 










r 


\bm{\it Gamma} 


A 


\bm{\itLaimbda} 


27 


\bm-[\ it Sigma} 




\bm{\itPsi} 


A 


\bm{\itDelta} 




\bm{\itXi} 


r 


\bm-T_\itUpsilon} 




\bm-[\ it Omega} 





\bm{\itTheta} 


// 


\bm{\itPi} 




\bm{\itPhi} 







Table 3: Boldface variants of slanted greek letters 



a 


\pmb-[ \upalpha} 


6 


\pmb-[\uptheta} 





\pmb{\upo} 


X 


\pmb-[\uptau} 


P 


\pmb-[\upbeta} 




\pmb-[\upvartheta} 


71 


\pmb{\uppi} 


t) 


\pmb-[ \upups i 1 on} 


Y 


\pmb-[ \upgamma} 


I 


\pmb-[\upiota} 


gj 


\pmb-[\upvarpi} 




\pmb-[\upphi} 


8 


\pmb-[\updelta} 


K 


\pmb{\upkappa} 


P 


\pmb-[\uprho} 




\pmb-[\upvarphi} 


£ 


\pmb{\upepsilon} 


X 


\pmb-[\uplambda} 


Q 


\pmb{\varrho} 


X 


\pmb{\upchi} 


e 


\pmb-[\varepsilon} 




\pmb-[\upmu} 


a 


\pmb-[\upsigma} 


V 


\pmb-[\uppsi} 


c 


\pmb{\upzeta} 


V 


\pmb-[\upnu} 




\pmb{\upvar sigma} 


0) 


\pmb{ \upomega} 




\pmb{\upeta} 




\pmb{\upxi} 










r 


\bm{\Gamma} 


A 


\bm{\Lambda} 


s 


\bm{\Sigma} 




\bm{\Psi} 


A 


\bm{\Delta} 


h- 1 


\bm{\Xi} 


T 


\bm{\Upsilon} 


ft 


\bm{\0mega} 





\bm{\Theta} 


n 


\bm{\Pi} 




\bm{\Phi} 







Table 4: Boldface variants of upright greek letters 



2 



